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Abstract. Very high energy gamma rays coming 
from extra-galactic sources can interact with inter- 
galactic radiation fields. This process may result 
in electromagnetic cascades with the following cy- 
cle: the production of electron-positron pairs and 
then secondary gamma-rays due to inverse Compton 
scattering. Since electrons and positrons will be 
scattered in the intergalactic magnetic field, under 
certain conditions their radiation may be redirected 
towards the observer. Thus one can anticipate that 
the secondary gamma-ray emission may produce an 
apparent extended halo around the source. 

MAGIC is an Imaging Atmospheric Cerenkov 
Telescope located on Canary island of La Palma at 
Roque de los Muchachos Observatory (2200 m.a.s.l). 
Various source sizes and extended emission profiles 
within 1° diameter have been studied by using 
dedicated Monte Carlo simulations for the MAGIC 
telescope. We present results of the study of a possible 
extended emission for Mrk 421 and MrkSOl done 
with the MAGIC telescope. 
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I. Introduction 

Blazars are well-known extragalactic emitters of Very 
High Energy (VHE) 7-rays. Their radiation transverses 
large distances in extragalactic space filled with CMB 
(Cosmic Microwave Background) and EBL (Extragalac- 
tic Background Light) photons and could be absorbed 
due to the pair production process. This effect leads to a 
change in the spectral shape of observed radiation and it 
can be used to constrain the EBL density (see e.g. [IJ). 

A possible extended emission from AGN objects was 
discussed in |2|. Gamma rays with energies > 100 TeV 
can be easily produced in hadronic models of AGNs. At 
those energies 7-rays will be strongly absorbed in a pair 
production process on CMB photons akeady relatively 
close to the source. For an Extra-galactic Magnetic 
Field (EGMF) strength of - 10 "^G, the produced 
e+e^ pairs will be almost immediately isotropized and 
efficiently cooled down in an inverse Compton process. 
This way a physical halo (of the radius > IMpc) of 
e+e^ pairs producing VHE 7-rays by Inverse Compton 
scattering can be created in the vicinity of the source. 
This extended emission is an input for the second step, 
where cascading on EBL photons takes place in the 
intergalactic space between the source and the observer 



Due to a relatively strong magnetic field, produced e+e^ 
pairs are scattered away from the observer field of view. 

According to most of the theoretical models a consid- 
erable fraction of the intergalactic space has to be filled 
with a very weak EGMF ([3|). 

Unlike Q the authors in Q, dU, Q are considering 
cascading in the much weaker EGMF of < lO^^^G. Be- 
cause of small deflection angles in weak magnetic fields, 
the e+e^ pairs created in the intergalactic space are 
not isotropized. In an inverse Compton process electron- 
positron pairs will produce secondary gammas with 
slightly different direction with respect to the primary 
photons. An electromagnetic cascade can develop if the 
optical depths for those secondary photons are still large 
enough . Secondary photons produced in those cascades 
can be deflected and redirected into the field of view of 
the observer. Also this process can lead to an apparent 
extended emission around the point-like source. 

The 7-ray emission from blazars is strongly variable. 
The variability timescales could be as short as a few 
minutes. In the case of an extended emission, depending 
on the distance from the source on which the pair 
production process occurred, the geometrical path-length 
transversed by 7-rays could be rather different. Due to 
this geometrical effect, the original short time variability 
will be smeared out to much longer time scales (^ 10^ 
years [4]). 

So far no extended emission around AGN object has 
been found (see e.g. [61). 

MAGIC is an Imaging Atmospheric Cerenkov Tele- 
scope located on Canary island of La Palma at Roque de 
los Muchachos Observatory (2200 m.a.s.l). Due to the 
large size of the mirror dish and improved light sensors, 
it has the lowest energy threshold among all lACTs. 
Recently the MAGIC project has been upgraded by 
adding a second telescope at 85 m distance for working 
in a coincidence (stereo) mode. 

II. Analysis method 

The image of an every given event can be 
parametrized by using the so-called Hillas parameters 
|7 |. The angular distance between the center of gravity 
of the image and the shower direction (so called DISP) 
is correlated with geometrical and timing properties of 
the image. Thanks to this the arrival direction for every 
event can be estimated. The distribution of the squared 
distance between the estimated and the true source 
position, (the 9^) is narrow and has a peak at = 
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Fig. 1. A compaiison of 9^ distributions obtained with parametrized 
DISP value (circles), and with a RF DISP (triangles) for the Crab 
Nebula data. 



for a point-like 7-ray source. In case of an extended 
source this distribution is broader 

For a single telescope the DISP is parameterized (see 
lis)). In a first approximation DISP is proportional to 
the ellipticity (1 — Width/ Length) of the image. By 
including the dependence on the image parameter Size 
and it's possible truncation because of the camera edge 
effects, one improves the precision of DISP. 

A. The novel Random Forest DISP method 

For this study we developed a novel method for the 
DISP estimation. As shown in [91 for a given image the 
time gradient along its main axis is strongly correlated 
with the impact parameter of the parent shower. Since 
for a given zenith angle and a 7-ray energy the DISP 
is a simple function of the impact parameter, one could 
expected that by using fast timing properties, one will 
improve the DISP estimation, and as a result also the 
angular resolution. 

To combine information from both geometrical and 
timing properties of the image in the most efficient 
way, we used multidimensional decision trees - the 
Random Forest (RF) method. It is widely used for 
the gamma/hadron separation and the energy estimation 
ITOl . A comparison of the novel RF method (hereafter 
RF DISP) and of the standard parametrized DISP is 
presented in fig.[T] The RF DISP provides a substantially 
narrower 9^ distribution. This improves the angular res- 
olution by ~ 20 — 30% and thus enhances the telescope 
performance for the search of an extended emission. 

The shape and the width of a 9^ distribution for a 
point-like source depends on many factors. The most 
important is the energy of the shower For higher energy 
showers, due to large number of particles in the shower 
maximum the relative internal fluctuations are smaller, 
and the resulting image has better defined parameters. 
Also the large Size (the total measured charge of the 
image) improves the precision of the reconstruction of 
the shower direction. In our analysis we used only 
showers with Size > 400photoelectrons, that provide 
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Fig. 2. Monte Carlo cumulative 8'^ distributions for a source with 
80% point-like and 20% extended emission. Characteristic radius of 
the extension is equal to 0.1° (dotted), 0.2° (dashed) or 0.3° (dot- 
dashed). The extended part of the emission is simulated as a flat disc 



(dN/de'^ 



const). Purely point-like source is shown with solid 



line. A random misspointing up to 0.03 deg has been included in 
the simulations. 



a relatively high precision of the 6^ determination. This 
leads to an energy threshold (defined as a peak of the 
Monte Carlo energy distribution) of 300 GeV. 

Let us consider a situation when the excess observed 
from a hypothetical source is a mixture of a point source 
and a small addition of extended emission with a given 
profile. The cumulative 9^ distributions for a point-like 
and extended sources are shown in fig. |2] Using this 
figure, the telescope PSF (point spread function) defined 
as 41% containment radius (equivalent to one standard 
deviation of a two dimensional gaussian distribution) 
of the reconstructed direction of 7-rays for the used 
analysis energy threshold is estimated to be < 0.1°. 

If the characteristic extension 9^x1 is smaller than the 
telescope's PSF, the distributions of a purely point-like 
and a partially extended source are very similar For 
larger extensions clear differences between distribution 
shapes can be seen (see fig. |2]i. 

To investigate possible extended emission from 
blazars we adopted a method similar to the one used 
in |6|. We calculated the ratio of event rates in two 9 
ranges; / = ^''(o^< /< ; ™d compared them with the 
ratio calculated for a point-like source. 

The 9i and 92 values are calculated with the help 
of Monte Carlo simulations searching for the most 
significant difference between purely point source and 
a source with an extended emission for a given size and 
profile of extended emission. 

To minimize the systematic errors we used the Crab 
Nebula as a model of a point-like source. This can be 
done because the mean spectral slopes of Mrk 501 and 
Mrk42I in the considered data sample are rather similar 
with that of the Crab Nebula. 

For every considered flux, radius and profile of the ex- 
tended emission by using MC simulations we calculated 
the value of / and of the corresponding significance of 
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Fig. 3. Difference between tlie true and the estimated source positions 
for Mrk421 (squares) and Crab Nebula (triangles). Circles with radii 
0.02° and 0.03° show the characteristic misspointing scale 



the extension. 

B. Misspointing 

One of the factors which limit lACTs ability to 
distinguish an extended from a point-like source is the 
misspointing of the telescope. Using a strong source like 
Crab or Mrk 421 one can estimate the source position 
and then calculate the misspointing as a difference 
between the true and the estimated positions. As it is 
shown in fig. |3] the characteristic misspointing scale is 
~ 0.02 — 0.03° This value is well below the angular 
resolution of the telescope and the assumed extensions 
of the halos. 

III. Data sample 

All data presented in this paper has been taken with 
a single telescope. To minimize systematic errors, only 
data taken at low (< 30°) zenith angle were used. Data 
from two different blazars have been analyzed. 

Data from Mrk 501 have been taken in April/May 
2008 in the so called ON/OFF mode (source is in the 
center of the camera). After quality cuts 26h of ON 
data were obtained. Estimation of the background is 
performed using 50h of OFF data. Mrk 501 in the above- 
mentioned time period was in a rather low state (mean 
flux in the sample 15% C.U. (Crab unit)) 

Mrk 421 data are from the period Dec 2007 - Jan 
2009. After quality cuts 32h of data were selected. 
Those data were taken in the so called wobble mode. 
Source position was shifted to 0.4° from the center of 
the camera. The opposite (with respect to the camera 
center) position was used for estimating the background. 
Since the background was taken simultaneously with 
the source data, the systematic errors were small. A 
disadvantage of this approach is the fact that if the 
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Fig. 4. The comparison of the 9'^ distribution for MrkSOl (circles) 
and a point-like source (triangles. Crab Nebula). 
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Fig. 5. The comparison of the 9^ distribution for Mrk421 (circles) 
and a point-like source (triangles. Crab Nebula). 



source extension is as large as 0.4°, the signal and 
the background regions start overlapping. This effect, 
that reduces the sensitivity was taken into account. In 
the analyzed data sample Mrk 421 was in a high state 
(- 1.4 C.U.). 

The 9^ distribution for a point like source was calcu- 
lated by using the Crab Nebula data. 

IV. Results 

0'^ distributions obtained for Mrk 501 and Mrk 421 are 
presented in the fig. |4] and fig. |5] For both sources the 
9^ distributions match with a corresponding point-like 
distribution. Using the first 12 bins, which contain most 
of the excess events we calculated /ndof=^OJ/ll (for 
Mrk 421) and 3.8/1 1 (for Mrk 501). 

An upper limit on the flux of the extended emission 
calculated for different extension radii and profiles is 
shown in fig |6] 

V. Conclusions 

We have developed a novel method based of the 
Random Forest for the estimation of the source posi- 
tion. This method improves the angular resolution by 
- 20 - 30%. 
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Fig. 6. An upper limit on tlie flux of tlie extended emission from Mrk 501 (left figure) and Mrk 421 (light figure) in the Crab units for 
different source profiles and extensions dN/dd (X0^: 0° < 9 < dcut- flat (/3 = 1) (open circles), 0.1° < d < Ocut- /3 = (stars), /3 = -1 
(crosses), f) = —2 (dots). 



For the used data samples stemming from 26h of 
observations of Mrk 501 and 32h of Mrk 421 no ex- 
tended emission has been detected around these sources. 
The study showed that if there is an extended emission 
around Mrk 501, then its flux is < 4% C.U. (see fig. |6] 
left), for the analysis energy threshold 300 GeV. In the 
case of Mrk 421, since the source was in a high emission 
state, the search for the halo was less sensitive, because 
point-like emission creates additional background. The 
constrain on the extended emission flux from Mrk 421 
is < 7% C.U. (see fig. fright). 

With the second telescope, the MAGIC system can 
operate in the stereo mode. By combining simultaneous 
information from both telescopes, the angular resolution 
will be significantly improved. This will allow us to per- 
form a search for the AGN halo with a better sensitivity. 

A future 7-ray telescope project - the Cerenkov Tele- 
scope Array is aiming for a much better angular reso- 
lution {^2 arc min). In this case lower halo extensions 
can be investigated with an improved sensitivity. 

Acknowledgements 

We would like to thank A. Neronov and D. Semikoz for 
drawing our attention to this interesting subject. 

References 

[1] Mazin, D., Raue, M. 2007, A&A, 471, 439M 

[2] Aharonian, E A., Coppi, RS., Volk H.J., 1993, ApJ, 423, L5 

[3] Elyiv, A., Neronov, A., Semikoz, D., arXiv:0903.3649 

[4] Neronov A., Semikoz D.V., 2007, JET? Lett. 85, 579 

[5] Dolag, K., Kachelriess, M., Ostapchenko, S., Tomas, R., 

arXiv:0903.2842 
[6] Aharonian, E A., et al. 2001, A&A, 366, 746 
[7] Hillas, A. M., 1985, Proc. 19th ICRC, La Jolla 3, 445 
[8] Domingo- Santamaria, E., Flix, J., Scalzotto, V., Wittek, W., Rico, 

J., 2005, Proc. 29th ICRC, Pune, 101 
[9] E. Aliu et al., Astropart. Phys. 2009, 30, 293 
[10] Albert, J. et al., 2008, NIM A, 588, 424 



